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Matching factors for AS = 1 four-quark operators in RI/SMOM schemes 
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The non-perturbative renormalization of four-quark operators plays a significant role in lat- 
tice studies of flavor physics. For this purpose, we define regularization-independent symmetric 
momentum-subtraction (RI/SMOM) schemes for AS = 1 flavor-changing four-quark operators and 
provide one-loop matching factors to the MS scheme in naive dimensional regularization. The mix- 
ing of two-quark operators is discussed in terms of two different classes of schemes. We provide a 
compact expression for the finite one-loop amplitudes which allows for a straightforward definition 
of further RI/SMOM schemes. 



I. INTRODUCTION 

The study of physical processes which change the 
strangeness by one unit (AS — 1), such as the decay of 
a kaon into two pions, is important for the understand- 
ing of CP violation within the Standard Model (SM) and 
its possible extensions. Such processes can be used to 
measure the parameter of direct CP violation e'/e, to 
study the AI =1/2 rule, and to calculate long-distance 
contributions to Kq — Kn mixing and the parameter of 
indirect CP violation e The resulting constraints 

for the Cabibbo-Kobayashi-Maskawa (CKM) matrix el- 
ements allow for a precise test of the SM. The weak in- 
teraction which mediates these processes with change in 
the strangeness can be described by local four-fermion 
operators at low energy scales, where the character of 
the vector boson interaction is essentially point like, see 
Refs. [1-G3 and Refs. [ll|-[l3| for reviews. Matrix ele- 
ments which describe, e.g., two-pion decays of kaons can 
then be computed with the help of lattice simulations. 

In order to perform the renormalization of relevant op- 
erators in the lattice computation one can adopt a renor- 
malization scheme which is independent of the regula- 
tor. Such a scheme can then be implemented in both 
non-perturbative lattice calculations and continuum per- 
turbation theory. This allows for a conversion of lat- 
tice results to the modified minimal subtraction (MS) 
scheme which is not directly applicable in lattice sim- 
ulations. In Ref. 14] the non-perturbative renormal- 
ization (NPR) technique and regularization independent 
(RI) momentum-subtraction schemes were defined for 
this purpose. 

In the context of light up, down, and strange quark- 
mass determinations quark bilinear operators need to be 
studied for the NPR procedure. The required match- 
ing factors which convert the quark masses and fields 
from the RI schemes to the MS scheme are known up to 
three-loop order [l4T - [l7j in perturbative Quantum Chro- 
modynamics (QCD). The renormalization constants in 
these regularization independent momentum-subtraction 
(RI/MOM) schemes are determined at an exceptional 
momentum point of the considered amplitude. In the 
case of the quark bilinear operators the exceptional mo- 



mentum configuration is distinguished by the fact that 
no momentum leaves the operator. However, a lattice 
simulation with an exceptional momentum configuration 
for the renormalization constants is more disposed to ef- 
fects of chiral symmetry breaking [l8[ . Furthermore, in 
the RI/MOM scheme unwanted infrared effects exist and 
the matching factors show a poor convergence behavior. 
For these reasons a non-exceptional momentum configu- 
ration was proposed in Ref. [18| and the framework and 
concepts of new RI/SMOM schemes with a symmetric 
subtraction point were worked out in Ref. [19). The sym- 
metric subtraction point is characterized by the fact that 
a momentum leaves the inserted operator. 

The matching factors for the conversion of quark 
masses from the schemes with a symmetric subtraction 
point to the MS scheme were computed to two-loop or- 
der [I9I - I22I ] by considering the amputated Green's func- 
tions with insertion of the scalar or pseudo-scalar oper- 
ator. These schemes exhibit a better infrared behavior, 
and also the coefficients of the perturbative expansion of 
the matching factors are smaller. Therefore their use led 
to a significant reduction of the systematic uncertainties 
in the light quark mass determinations following this ap- 
proach [23| compared to previous studies, see Ref. [24| . 

For the insertion of any multi-quark operator into an 
amputated Green's function the fermion field for each 
external leg needs to be renormalized. In Ref. [l9| two 
schemes, the RI/SMOM and RI/SMOM^ scheme, were 
suggested for the renormalization with a symmetric sub- 
traction point. It was also shown that the former is 
equivalent to the RI'/MOM scheme and is thus known to 
three- loop order JT3 , whereas the latter is known 

to two-loop order 



1< 



- 22] . The corresponding calculation 
requires the computation of amputated Green's functions 
with insertion of the vector or axial- vector operator and 
the utilization of Ward-Takahashi-identities. 

We would like to mention that also the tensor op- 
erator has applications in lattice simulations, see, e.g., 
Ref. [25[ and a scheme with a symmetric subtraction 
point has been introduced in Ref. 19]. Its matching 
factor and anomalous dimension is known to two-loop 
order [l9|, [2l|, [22|. Also moments of twist- 2 operators 
used in deep inelastic scattering have been studied in a 
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RI/SMOM scheme in Refs. [2^, H3- 

In view of these successes and advantages, the 
RI/SMOM definition has been extended to AS = 2 
flavor-changing four-quark operators in Ref. (28[, where 
the one-loop QCD corrections to different matching fac- 
tors have been computed, and also the anomalous dimen- 
sions were provided. These results were then used for 
the determination of the Bk parameter which is needed 
to parametrize the hadronic matrix element for the the- 
oretical description of K — K mixing. For the case 
of an exceptional subtraction point the matching factors 
were determined at next-to-leading order in Refs. [291.130} 
based on Ref. 3l|. Similarly the matching for AS = 1 
flavor-changing four-quark operators with an exceptional 
subtraction point was determined in Refs. 0,|3|. The 
purpose of this paper is to introduce RI/SMOM schemes 
with a non-exceptional subtraction point for AS = 1 
flavor-changing four-quark operators as well as to pro- 
vide the corresponding matching factors for the conver- 
sion from these RI/SMOM schemes to the MS scheme in 
naive dimensional regularization (NDR). To this end we 
first present the framework needed to properly take into 
account the mixing with two-quark operators which was 
not needed in the AS = 2 case. We then study the in- 
sertion of the AS = 1 operators into amputated Green's 
functions in perturbative QCD at one-loop order to de- 
termine the renormalization constants. 

The outline of this work is as follows. In Sec. |TT] we 
define the set of AS = 1 four-quark operators used in 
this work. In Sec. IHII we discuss some generalities of the 
renormalization of the AS = 1 operators in the MS [NDR] 
scheme as well as in a general RI scheme and introduce 
our notation. In Sec. IIVI we provide a classification of 
projectors used to define the RI schemes and present our 
results for the finite one-loop amplitude as well as con- 
version factors from different RI/SMOM schemes to the 
MS [NDR] scheme. Finally we close with a summary and 
conclusions in Sec. El 



II. THE BASES OF AS = 1 OPERATORS 

In this section we define operator bases of the effective 
AS* = 1 Hamiltonian of electroweak interactions, where 
we closely follow the notation of Ref. [HI . We work in an 
effective three-flavor theory including the up, down, and 
strange quark. This effective theory is valid for energies 
below the charm quark mass. The effective Hamiltonian 
reads 



njAS=l 



(i) 



with Fermi coupling constant Gf and renormalization 
scale \i. The symbols Cj(/i) denote Wilson coefficients 
and Oi(fi) are four-quark operators, which we will discuss 
in terms of "physical" and "chiral" operator bases in dif- 
ferent schemes labeled x. In the case of the physical oper- 
ator bases the effective Hamiltonian is expressed in terms 



of ten operators which are grouped into current-current 
operators, QCD penguin operators, and electroweak pen- 
guin operators. The physical operator bases are classi- 
fied by the physical origin of their respective operators, 
whereas the chiral operator basis is classified by irre- 
ducible representations of the SU(3)l <8> SU(3).r chiral 
symmetry. 



A. The physical bases 

Let us start with the traditional physical operator basis 
of Refs. 13414361 ]. The current-current operators are 
defined by 



Qi = {s a Ub)v-A{ubd a )v-A 
Q-2 = (s a u a )v~A(ubdb)v-A 



(2) 



the QCD penguin operators are defined by 

Q3 = {s a d a )v-A ^2 {qbqb)v~A , 

q— u,d,s 

Qi = {s a db)v-A ^2 (QbQa)V-A , 
q— u,d,s 

Q5 = {s a da)V-A ^ (q~b<lb)v+A , 
q—u,d,s 

Qa = {s a db)v-A ^2 (Qb(la)v+A , (3) 

g— u : d,s 

the electroweak penguin operators are defined by 

Q7 = ^(Sad a )v-A 2J e q(%qb)V+A, 
q—u,d,s 

Qs = ^{s a db)v-A ^2 e q(^bq a )v+A (4) 

q— n,d,s 

with e u — 2/3, ea = e s = —1/3, and 



Q9 = ^( S ada)V-A ^2 e l(lbqb) 



V-A 



q—u,d,s 



Qw = -^{s a db)v-A ^2 e q(qbq a )v-A 

q—u,d,s 



(5) 



where {qq)v±A refers to the spinor structure <?7 M (l±75)g, 
a and b are color indices, and u, d, and s are the fields of 
the up, down, and strange quarks. This basis of operators 
{Oi} = {Qi, • • • , Q10} is referred to as "basis I" in the 
following. Alternatively we can Fierz transform Qi and 
Q 2 to 



Ql = {s a d a )v-A{ubUb)v-A 
Q 2 = {s a db)v-A{ubU a ) V -A 



(6) 



The basis of operators {0 L } = {Qi, Q2, Q3, ■ ■ ■ , Qw} is 
called "basis II" in the following. Operators with color 
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contractions as in Qi are called color diagonal, operators 
with color contractions as in Q 2 are called color mixed. 
It will also be useful to define the Fierz transformation 
of Q 3 , i.e., 

Q3 = ^2 (s a qb)v-A(qbd a )v-A ■ (7) 

q— u,d,s 

In an explicit four-dimensional regularization scheme, 
such as lattice regularization, Qi and Qi can be used 
interchangeably. In dimensional regularization, however, 
the contribution of evanescent operators such as 

E\i = Qi~ Qi (8) 

has to be included. The evanescent operators vanish at 
tree level if the regulator is removed, see, e.g., Refs. |35L 

MM- 

B. The chiral basis 



(8,8) Q' 78 = Q 7)8 , (10) 

where {L, R) denotes the respective irreducible represen- 
tation of SU(3)i<g)SU(3)fl. Instead of expressing Eqs. © 
and (fit))) in terms of the operators of basis I and evanes- 
cent operators we could have also eliminated the latter 
by introducing the operators Q\, Q2 of basis II and Q3 
with the help of Eq. ©. 

III. RENORMALIZATION 

In this section we discuss the renormalization of the 
four-quark operators Oi starting with some generalities 
concerning operator renormalization at fixed gauge. We 
then describe the one-loop off-shell renormalization of the 
operators in the different bases in the MS[NDR] scheme 
in detail and provide a discussion of renormalization in 
a general RI scheme. A brief discussion of some details 
concerning the Wilson coefficients of the chiral basis con- 
cludes this section. 



The operators Qi,...,Qio are not linearly indepen- 
dent, i.e., one can eliminate three operators by express- 
ing them as linear combinations of the remaining ones. 
In a regularization which breaks Fierz transformations 
also evanescent operators enter these relations. The re- 
duced operator basis of linearly independent operators 
can then be classified according to irreducible represen- 
tations of SU(3) and SU(2) flavor symmetries @, Hj], and 
the resulting operator basis will be referred to as the 
"chiral basis" in the following. The linear independence 
of its elements will become important later for the non- 
perturbative definition of RI schemes with the help of 
projectors. 

We proceed along the lines of Ref. [2J amending their 
discussion by the contributions of evanescent operators 
since we work in dimensional regularization. We first 
eliminate the operators Q4, Qg, and Qiq using 

Q4 = Q2 + Q3 — Qi — E12 — Ei 3 , 

3 13 
Q9 = 2<9i _ _ 2 En ' 

Qw = \{Qx-Q3) + Q2 + \e^-E 12 . (9) 

The remaining seven operators can then be recom- 
bined according to irreducible representations of the chi- 
ral flavor-symmetry group SU(3)l<S>SU(3)_r. The details 
of this decomposition including evanescent operators are 
given in App. [S] The chiral operator basis is thus given 
by 



~(2Qi- 2Q 2 + Q 3 ) - g£i 
i(-3Qi + 3Q 2 + Q 3 ) + ^ 

5,6 = Qbfi , 



(27,1) 


Q'i 


(8,1) 


Q'i 


(8,1) 


Qs 


(8,1) 


Q' 5fi 



A. Renormalization at fixed gauge 

The renormalization schemes described in this work 
are defined at a fixed covariant gauge with gauge fix- 
ing parameter £, where £ = (£ = 1) corresponds to 
the Landau (Feynman) gauge. The gauge-fixing proce- 
dure explicitly breaks the gauge symmetry, and it can 
be shown that mixing with three classes of operators can 
occur [35l . |40M44| : (i) gauge-invariant operators which 
do not vanish using the equations of motion, (ii) gauge- 
invariant operators which vanish using the equations of 
motion, and (iii) gauge non-invariant operators which are 
either BRST invariant or vanish using the equations of 
motion. 

In a general renormalization scheme x, a renormalized 
four-quark operator Of can be written as 

of = y.;p, + b* k F k + c^d + di n N m , (ii) 

where Oj are the bare four-quark operators, F k are 
evanescent operators, Gi (N m ) are gauge-invariant 
(gauge non-invariant) operators involving only two quark 
fields. The symbols Zfp bf k , cf ; , and df m denote renor- 
malization constants. The sum over the respective op- 
erator basis for Oj, F k , G\, and N m is implied. The 
operators Oj belong to class (i), the operators Gi belong 
to class (i) or (ii), and the operators N m belong to class 
(iii). Operators of clas s (ii) a nd (iii) do not contribute 
to physical amplitudes (40l - |43l | . The mixing of operators 
Nm can be avoided by using the background field gauge 

fallal HI. 



B. The MSJNDR] scheme 

In the following we discuss the off-shell renormaliza- 
tion of the operators Oi in the MS scheme with massless 



quarks at one-loop order in perturbative QCD. We use 
naive dimensional regularization (NDR) in d = 4 — 2e 
space-time dimensions with a naive anti-commutation 
definition of 75. For multi-loop calculations the opera- 
tor basis of Ref. [46j is more convenient and allows for 
a straightforward treatment of 75. Since we restrict our- 
selves in this work to the one-loop order, we adhere to 
the traditional bases in order to connect with previous 
works in this field. The MS-renormalization of the four- 
quark operators Qi with a focus on on-shell renormaliza- 
tion is discussed, e.g., in Refs. [Hil!?} at one loop and in 
Refs. [H, HI] at two-loop order. 

The four-quark operators Oi are of mass-dimension 6, 
so that in the massless limit mixing can only occur with 
operators Fk, Gi, and N m of mass-dimension 6. For off- 
shell external states the operator 



Gt 



;S lv {l- lb )[D^[D^D»]]d 



(12) 



mixes under renormalization with the operators Oi at 
the one-loop level [H, HU . This operator is of class (i) , 
i.e., it is nonzero in the limit of on-shell external states. 
In this limit, however, the operator G\ becomes linearly 
dependent on the four-quark operators Oi , and one finds 

Emm 



Qp 



with 



(13) 



(14) 



At one-loop order the renormalized operator in the MS 
scheme is given by 



, MS , 



5 E 1 _1_ J4S^ 



with 



MS 



zMS 



(15) 



(16) 



There are two different types of diagrams that need 
to be considered: current-current diagrams, where each 
fermion line involving a quark field of the operator Oi 
extends to the external quark fields, and penguin dia- 
grams, where one fermion line involving quark fields of 
the operator Oi begins and ends at the operator. We de- 
pict the corresponding one-loop diagrams in Fig. [1] and 
Fig. [5J The one-loop current-current diagrams determine 
the mixing of the four-quark operators with themselves 
(given in a MS ), while the penguin diagrams determine 
the mixing of G\ with the four-quark operators (given in 

We first consider the off-shell renormalization of oper- 
ator basis I and II. The current-current contributions in 



,MS 



separate in 2 x 2 blocks given by [10|, |3J, |36|, |47 




FIG. 1. Current-current contributions at one-loop order. 




FIG. 2. Penguin- type diagrams with up to three external 
gluons. Three analog penguin diagrams, where the gluons 
attach to a closed fermion loop, also need to be considered 
but are not shown here. 



47T6 



-3/N c 3 
3 ~3/N c 



and 



MS -1 



47T6 I 



(3/N c 

3/N c 



3N t 



(17) 



(18) 



where the subscripts i, j of (cqj ) denote that the matrix 
acts on the space of operators Oi and Oj . This notation 
is also used for other block-diagonal matrices in the re- 
mainder of this section. The matrix a MS is identical for 

basis I and II. 

The penguin contributions in c MS are given by 



1 



5 



JVIS _ „MS 



47T£ 



.MS 

'1,5,7,8,10 



(19) 



for basis I and II. For on-shell matrix elements we use 
Eq. (Tl"5|) in Eq. (fT5|) and thus reproduce the results for 
the anomalous dimensions of Refs. [ltj HH [36ll47|. 

The set of evanescent operators {Fk} used in Eq. (fT5)) 
to define the MS scheme consists of operators 



Ei = {-Sal^l L p u b ){u b y i YY L d a ) - (16 - 4e)Qi 

El = {Salnlulp U a ){u b Y-f V Y L d b ) 



- (16-4e)Q a . 



(20) 



and 



Ei = {s allll , 1 L p d a ){u b Y 1 v Y L u b ) - (16 - 4e)Qi , 

El = {Salplulpd b )(u b YYl PL Ua) 

-(16-4e)Q 2l (21) 
as well as 

q=u,d,s 

- (16 - 4e)Q 3 , 

Ei = {sainivipd b ) ^2 (qbj tl j l/ i pL q a ) 

q—u,d,s 

-(16-4e)Q 4 , 
£ 8 = (wvtfda) E fe7 M 7V%) 

q— u,d,s 

- (4 + 4e)Q 5 , 



q—u,d,s 



(4 + 4e)Q 6 



(22) 



and 



E 7 = -{salplvl L p d a ) ]T e 9 (& 7 'W B 9&) 



q—u,d,s 



(4 + 4e)Q 7 , 



E & = -{s al p lvl L p d b ) Y, zMl P l y l pR 1a) 

q—u^d^s 

-(4 + 4e)Q 8 , (23) 



and finally 



^9 = y(Salplvlpd a ) Y e q(qbYYl PL 1b) 



q—u,d,s 



(16 - 4e)Q 9 



#10 = ^(Salp7u7pd b ) Y e q {q b Yl V l pL qa) 
q—u,d,s 

-(16-4e)Qi 0) (24) 



where 7p ^ = 7 P (1 T 7s), an d a and 6 are color indices. 
The explicit contributions of order s are determined by 



the "Greek" method 14911 in accordance with two-loop cal- 
culations such as Ref. [30]. For operator basis I we have 
{Fk} — {Ei, . . . ,Eiq}. The renormalization coefficients 
of the evanescent operators in Eq. (|15[) decompose in 2 x 2 
blocks and are given by 



(bfl) = 



1/4 
-l/(2iV c ) 



(N e /4-l/(2N, 
Aire \ 1/2 

(bW = (bW) - = (fiffo) 

as (-\/{2N c ) 1/2 
4tt £ I 1/4 iV c /4 - 1/(2JV C ) 



(25) 



for basis I. Similarly for basis II we have {F^} = 
{Ei, Ei2, Es, ... , Eiq}, and one obtains 



(b 



MS\ 
3,4 I 



4tt£ 



= («#f) = Kl) = «1 S 0) 

-l/(2iV c ) 1/2 
1/4 iV c /4-l/(2JV ) 



(26) 



If operators transform in an irreducible representa- 
tion of a given symmetry, they only mix with other op- 
erators transforming in the same irreducible represen- 
tation. In the case of the chiral basis the decompo- 
sition of operators according to irreducible representa- 
tions of SU(3) L <g> SU(3) fl is given in Eqs. JT0}. The 
set of operators used in Eq. ([To]) for the chiral basis is 
given by {O t } = {Q[, Q> 2 , Q' 3 , Q> 5 , . . . , Q> 8 } and {F k } = 
{Ei, . . . , E13}, where the operators En, E12, and £13 
are defined in Eq. (JSJ. The matrix a MS is then again 
block-diagonal with 

(af S ) = ^(3-3/N c ), 

r„MS\_ a ° (~Z/ N c 3 
[a2 > 3) ~ Aire { 3 -3/N c 

K, 6 ) - K,s) - AjT£ y o 3/jv c -3JV c '' [ -" 



and 



47T£ V 3 ' 



(28) 



The nonzero elements of b are given in Tab. UJ In the 
chiral basis the on-shell limit of Gi is given by 

Q p = (2 - 3/N c )Q' 2 + (3 - 2/N c )Q' 3 - Q' 5 /N c + Q' 6 

— En — E12 — Ei 3 . (29) 



C. Regularization-independent schemes - The 
mixing of four-quark operators 

In the following we define RI schemes for the AS = 1 
four-quark operator bases. The RI schemes are defined 
non-perturbatively, so that they can be used in lattice 



6 





jMS / Q s 




L MS / a s 


(1,1) 


1 


(3,1) 


3/10 


(1,2) 


-1/N C 


(3,2) 


-3/(10iV c ) 


(1,9) 


-1/JVc 


(3,9) 


1/(5JV C ) 


(1,10) 


1 


(3, 10) 


-1/5 


(1,11) 


9 


(3,11) 


6/5 


(1,12) 


-6 


(3, 12) 


6/5 


(1,13) 


3 


(3, 13) 


-3/5 


(2,1) 


-1/5 


(5,5) 


-l/(2JVc) 


(2,2) 


1/(5JV C ) 


(5,6) 


1/2 


(2,3) 


-1/(6AT C ) 


(6,5) 


1/4 


(2,4) 


1/6 


(6,6) 


N c /A - 1/(2JV C ) 


(2,9) 


-2/(15iV e ) 


(7,7) 


-V(2JVc) 


(2, 10) 


2/15 


(7,8) 


1/2 


(2,H) 


-9/5 


(8,7) 


1/4 


(2, 12) 


-9/5 


(8,8) 


iV c /4 - 1/(2JV C ) 


(2, 13) 


-3/5 







TABLE I. The nonzero elements of 6 for the chiral basis. 



simulations as well as in continuum perturbation theory. 
In lattice calculations the RI schemes serve as intermedi- 
ate schemes and allow for a straightforward conversion of 
the studied quantity to the MS scheme. In this subsec- 
tion we focus on the mixing of four-quark operators Oi 
among themselves. In terms of Eq. this means that 
we provide the RI renormalization conditions to deter- 
mine the renormalization matrix Z Rl . The renormaliza- 
tion conditions which determine the mixing of two-quark 
operators Gi and N m with the four-quark operators Oi, 
i.e., c RI and d m , are discussed in subsection MI Dl While 
the content of this subsection thus suffices to define the 
RI schemes for the (27, 1) and (8, 8) operators of the chi- 
ral basis, the discussion of subsection MIDI is necessary 
to complete the RI schemes for the (8, 1) operators. 

Let us consider a set of bare operators {Oi} that is 
closed under renormalization and contains the set of four- 
quark operators {Oi}, i.e., {Oi} C {Oi}. The renormal- 
ized operators Of in the MS or an RI scheme can be 
expressed in terms of the bare operators by 



of z;,o,. r,. {ms,ri}. 



(30) 



The renormalization conditions of the RI schemes are 
formulated in terms of renormalized amputated Green's 
functions r^(Of) with a single insertion of such an op- 
erator Of, where y denotes the wave function renormal- 
ization scheme and n indicates the external states of the 
Green's function. In this subsection we only consider 
Green's functions with four external quarks which we de- 
note by n — 4. 

The renormalized amputated Green's function T\ is 
related to the bare amputated Green's function T4 by 



ri(oi 



i 



(z y q ) 



(31) 



where Z y is the quark wave function renormalization con- 
stant, which relates the bare quark field / to the renor- 

1/2 

malized quark field p = (Z y ) 1 f in the scheme y. 

Various RI/SMOM wave function renormalization 
schemes have been proposed in Ref. [l9[ and will be used 
later. In order to convert the quark fields from the RI 
scheme to the MS scheme matching factors 



yMS 

A/ 



(32) 



have been computed with / MS = {C y q ) 1/2 p . 

The operators which are renormalized in an RI scheme 
can be converted to the MS scheme 



O 



MS 



gRI-s-MS^jRI 



using the conversion matrix 



S 



RI->MS 



zMS/^RI 



(33) 



(34) 



The renormalization matrix Z has been discussed in 
the previous Sec. IIIIB1 whereas the matrix Z m is de- 
termined by the renormalization conditions of the RI 
scheme. Using Eqs. (|3"2"]l and (|3~4"|) one can express 
Eq. (|31j) completely in terms of matching factors and 
renormalized Green's functions 



rl(of) 



RI^MS 



Hi L 4 [ U i 



(35) 



In our case the set of operators {Oi} is given by 
{Oi,Gi,N m ,Fk}, where {Oi} are the four-quark opera- 
tors of basis I, II, or the chiral basis and {-Ffc} are the 
corresponding evanescent operators used to define the 
MS scheme in the previous section. Since the evanescent 
operators are an artifact of dimensional regularization 
and escape a regularization-independent definition, their 
contribution to the right-hand side of Eq. (|3"3")) should be 
avoided. In order to achieve this, one defines new sub- 
tracted operators {Of ub } from the set of bare operators 
{O i } = {O i ,G l ,N m }. 

In general for any regulator one first has to subtract all 
contributions specific to the regularization from a given 
bare operator. Hence, in dimensional regularization we 
perform modified minimal subtraction of the evanescent 
operators [3^,13^, i.e., the subtracted operator is defined 
as 



0, sub = Oi 



.MS 



(36) 



where s ° is chosen such that it cancels all contributions 
of Fk proportional to a pole in e. In principle the choice 
of evanescent operators {Fk} used on the right-hand side 
of Eq. (|36|) is not unique. A useful choice is the basis {Fk} 
given in the previous section to define the MS scheme for 
operator basis I, II, and the chiral basis, and for which 
we therefore have 



O 



sub, MS 



O 



MS 



(37) 
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For convenience we adopt this definition of the subtracted 
operator in the following. In a lattice regularization one 
has to perform a similar subtraction of lower-dimensional 
two-quark operators Q which do not occur in dimen- 
sional regularization. From now on we only consider such 
subtracted operators and therefore drop the explicit no- 
tation of the superscript "sub" . 

In the RI schemes one imposes the renormalization 
condition |14j that amputated Green's functions with 
given off-shell external states at a given momentum point 
and in a fixed gauge coincide with their tree-level value. 
This condition is made explicit by choosing a certain set 
of projectors {Pij} in spinor, color, and flavor space that 
is applied to the four-quark amputated Green's function 
and imposing 

p ^liof l )\ m conf = p 4j rr c (o.) , (38) 

where r 4 roo (Oi) denotes the insertion of the operator 
Oi in the amputated Green's function r 4 at tree level. 
Since we are only interested in Rl-renormalized opera- 
tors Of 1 , we do not provide RI conditions for the op- 
erators Gi and N m in this work. For a set of n oper- 
ators Oi G {0%, . . . ,O n } we will provide n projectors 
Pij G {Pa, . . . , Pin] to determine the n x n elements of 
the renormalization matrix 2T RI . If no two-quark opera- 
tors mix with the four-quark operators Oi, i.e., c RI = 
and d RI = in Eq. the renormalization matrix is 

given by 

Z m = (Z^) 2 F 4 [M 4 (0)]- 1 , (39) 

where 

[Mi(0)] VJ = P 4j T 4 (O l )| mom conf • (40) 

and Fi = lim Qs ^o M\{0). If two-quark operators mix 
with the operators Oi, Eq. (I3T))) has to be modified 
slightly, which is discussed in subsection IIII D[ 

The matrix Z m given in Eq. (|39| depends on the reg- 
ulator used to define the RI scheme, i.e., the matrix Z RI 
obtained using a lattice regulator is different from the 
matrix Z RI obtained in dimensional regularization. The 
RI condition of Eq. (|38j) . however, fixes the physical am- 
plitudes of the Rl-renormalized operators at a certain 
off-shell momentum point to its tree-level value, which 
is independent of the choice of the regulator. Therefore, 
the physical amplitudes of the Rl-renormalized operators 
agree for all choices of the regulator. 

We define the RI scheme in the limit of vanishing quark 
masses. The choice of projectors {Pij}, the gauge fixing, 
and the momentum configuration of the off-shell ampu- 
tated Green's functions T\ defines the scheme up to mix- 
ing with two-quark operators. The explicit form of the 
projectors will be discussed later in Sec. IIVI Note that 
Eq. ([55)1 matches the amputated Green's function of a 
certain physical process with insertion of RI operators at 
a certain off-shell momentum point. 

In particular for the AS = I operators we consider the 
off-shell amputated Green's functions Y^ 3 ^ 8 ' 1 ^ 1 ^ of the 



ff(Pl) IU-P2) 




2q 



4V) S T(-^) 

FIG. 3. Momentum flow and open indices for the four-quark 
Green's functions r 4 . We explicitly write the quark fields sf , 
dj, f^, and }f with spinor indices a, /3, 7, and 6 and color in- 
dices i, j, k, and I. The momenta of the quark fields are given 
in brackets, which are counted as incoming, i.e., they flow to- 
wards the four-quark operator. An additional momentum of 
2q leaves the operator as indicated by the arrow. 



process 

d( Pl )s(- P2 ) -> f(- P3 )f(pi) (41) 

with quarks f = u,d, s, momenta p%, p 2 , P3, Pi, color in- 
dices k, I, and spinor indices a, /3, 7, 8. Using crossing 
symmetry we could equally well consider the scattering 
amplitude 

d(pi)/(ps) "> «(P2)/(P4) ■ (42) 

The momentum configuration used in Eq. (1381) to define 
the RI/SMOM scheme is then given by 

P3 = Pi , Pi= P2 (43) 

with 

222 2 I A A\ 

Pi = P2 = 1 = "Ms i Q = Pi~P2, (44) 

in Minkowski space, where /i s is the subtraction scale. 
This momentum configuration and our convention for 
open indices for the amputated Green's functions T 4 is 
shown in Fig. [3] In Fig. @] we explicitly show the cor- 
responding penguin diagrams at one-loop order, where a 
momentum transfer 2q leaves the operator. This choice of 
momenta is non-exceptional (no partial sum of incoming 
external momenta vanishes) which has the advantage of 
suppressing unwanted infrared effects in the lattice sim- 
ulation. 

In Ref. [13] a RI/MOM scheme was defined which uses 
exceptional kinematics and a different momentum point 
for current-current and penguin diagrams, see Fig. 5 of 
Ref. [32| . The momentum configuration in the RI/MOM 
scheme is given by p — pi — pi = P3 = Pi at fij. = —p 2 
for current-current diagrams and p = p\ — Pi and 
p' = P2 = P3 at fi 2 . = —q 2 for penguin diagrams. We 
consider the RI/MOM scheme in the following for com- 
pleteness, illustration, and check of our calculation. An- 
other RI scheme with exceptional kinematics which uses 
the same momentum point for current-current and pen- 
guin diagrams is discussed in Ref. (Hoj . 
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ffiPl) JU-P2) 

* \ * 

O 

o 
o 
o 
o 

^^■u,d,s 
> O > 

d?(pi) 

(6) 

FIG. 4. Penguin contributions at one loop and the respective 
momentum configuration. We explicitly write the fields and 
momentum configuration as in Fig. [3] We do not explicitly 
show the flow of the additional momentum 2q here. 



D. Regularization-independent schemes - The 
mixing of two-quark operators 

In the following we discuss the mixing of the two-quark 
operators Gi and N m with the four-quark operators Oi in 
the RI schemes. Such mixing occurs, e.g., for the (8, 1) 
operators of the chiral basis. The two-quark operators 
mix through the penguin diagrams, and therefore their 
mixing should be determined from amplitudes which only 
receive contributions from penguin-type contractions. 

Two kinds of such amplitudes will be considered in 
the following: (a) amputated Green's functions T 2 with 
two external quarks and one external gluon and (b) am- 
putated Green's functions T^ p with four external quarks 
corresponding to the process df — > sf, where the quark 
flavor / ^ {u, d, s}. The momentum flow and our conven- 
tion for open indices in case (a) is shown in Fig. [5] The 
corresponding momenta px, P2, and q for the RI/SMOM 
scheme satisfy Eq. (|44p . This momentum configuration 
is also non-exceptional. In case (b) we adhere to the 
momentum configuration and choice of indices shown in 
Figs. [3] and 2J In the following we define the RI schemes 
in such a way that both cases, (a) and (b), lead to iden- 
tical one-loop conversion factors from the RI scheme to 
the MS[NDR] scheme. At higher loops, however, the RI 
schemes defined by (a) differ from the RI schemes defined 
by (b). Case (b) can be implemented in lattice simula- 
tions in a straightforward way Q , while case (a) requires 



FIG. 5. Momentum flow and open indices for the two-quark 
Green's functions F2. We explicitly write the quark fields 
sf, dj and the gluon field A^ with spinor indices a, /?, color 
indices i, j, and Lorentz index fi. The index a = 1, . . . , N% — 1 
is contracted with the generators of SU(iV c ). The incoming 
momenta of the fields are also given in brackets. 

an external gluonic state. Nevertheless, the availability 
of both cases should be beneficial for lattice simulations, 
in particular to estimate higher-loop effects that are ne- 
glected in this work. 

Case (a) is formulated in terms of renormalized ampu- 
tated Green's functions T| which are related to the bare 
amputated Green's functions T 2 by 

n?(Of) = zg( ^ )1/a g5r a (o J ), (45) 

where Z y A is the gluon wave function renormalization con- 
stant, which relates the bare gluon field to the renor- 
malized gluon field A^ = {Z v A ) 1 / 2 A tl in the scheme y. 
Using Eqs. (|32|) and ([34| one can express Eq. (|45l> as 

r|(0f) = (zT/zi^cy 

x (5 RI ^)^. 1 r^(Op) . (46) 

We then determine the renormalization coefficients c RI 
and d RI by imposing 

p 2fc r|(of I )| mom conf . = p 2fe r t 2 rcc (o J ) = (47) 

for a certain set of projectors {P 2 k}- Since only tree- 
level insertions of Gi and N m need to be considered in 
this work, we do not provide renormalization conditions 
to define Rl-renormalized operators Gf 1 and iV RI . The 
definition of RI schemes for the two-quark operators Gi 
and N m is beyond the scope of this work. 

The RI conditions given in Eqs. (J38]) and (g7|) then 
allow for a non-perturbative determination of the renor- 
malization matrices Z RI , c RI , and d m defined in Eq. 
Without loss of generality we set d RI = in the following. 
We find 

Z m = (^) 2 F 4 [M 4 (0)]- 1 (1-^ RI )- 1 , 
c RI = -Z m M 2 {0)[M 2 {G)]- 1 (48) 
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with 



and 



M 2 (O) [M 2 (G)] _1 M 4 (G) [M 4 (O)} 



[M 2 {0)] ik = P2kT 2 {O l 



(49) 



(50) 



Case (b) is formulated in terms of the renormalized am- 
putated Green's function r| p corresponding to the pro- 
cess sf — > df that is shown in Fig. 2] at one loop. Note 
that the quark flavor / {u, d, s} for T v ip as opposed to 
/ € {u,d, s} for T V A which was used in subsection IIII CI 
We then determine the mixing of the two-quark operators 
by imposing 



= Pi P ,kVX e {Oi) = (51) 



for a certain set of projectors {-P4 P) fc}. The renormal- 
ization coefficients for case (b) can be obtained from 
Eqs. (J35J) and fU]) by replacing M 2 



M^ p with 



[M 4p (0)] ik = P 4p!k T 4p (Oi)\ 



(52) 



We provide the explicit projectors used to define the 
RI/SMOM schemes in Sec. [[V] We will define the set 
of projectors {P 2 k} and {Pi p .k\ such that the resulting 
RI/SMOM schemes agree at one loop. From now on 
we refer to RI/SMOM schemes defined using case (a) as 
RI/SMOM 2 schemes, while we do not write a subscript 
for RI/SMOM schemes defined using case (b). 

Since in the physical application of the RI scheme one 
is not interested in the conversion of the operators Gi and 
N m from the RI scheme to the MS scheme, it is useful to 
decompose the matrix S^ 1- *' s of Eq. ([34]) into several 
blocks. The conversion relation of Eq. ([55)1 for subtracted 
operators is then given by 



a MS = o 



RI 



Ad 



L J J 

RI^MS^yRI 



Ac RI^MS G RI 



(53) 



where the operators Oi are either of basis I, II, or the 
chiral basis. The three blocks Aa^ MS , Ac RI ^ MS , and 

Adf^ MS are obtained by imposing the renormalization 
condition of Eqs. (j3"8"j) and (|4"7|) . i.e., we have to evaluate 
r Ms( p) with of j. slie u external legs, see Eq. p]) . 

We denote on-shell matrix elements of the opera- 
tors Of 1 with a general external momentum setting by 
Oj^ s > for which one finds 



O 



MS\ 



R 



RI- 



■MS ^QRl^ 



(54) 



since the on-shell matrix elements of the two-quark op- 
erators Gi and N m are related to the on-shell matrix ele- 
ments of Oi . At the one-loop level only the two-quark op- 
erator G\ contributes to the right-hand side of Eq. ((52 
and one has 



(Gi) = (Q P ) ■ 



(55) 



In Sec.|TV]we compute Aa RI ^ MS and Ac RI ^ MS as well 
as the on-shell conversion factors 



Ar 



RI^MS 



R 



RI-yMS 



(56) 



at one loop for different RI schemes and for operator basis 
I, II, and the chiral basis. 



E. The Wilson coefficients of the chiral basis 

For lattice calculations it is advantageous to consider 
the chiral operator basis that uses the classification of 
operators according to irreducible representations of the 
chiral symmetry. In this basis the operators of each ir- 
reducible representation can be renormalized indepen- 
dently. Therefore the effective Hamiltonian should be 
expressed in terms of the operators Q'f I . 

On-shell matrix elements of the effective Hamiltonian 



defined in Eq. (JT|) in terms of Q \ read 



= ^J2Cf%)R' 



RI-i-MS 
ij 



(m) 



(57) 



where /i is the renormalization scale and n ,Kl ~^ MS j s 
conversion matrix of Eq. (|54[) for the chiral basis. The 
Wilson coefficients are, however, typically given for the 
traditional operator bases in Refs. (33,HiI[5l[[5|j]- There- 
fore it remains to relate the above Wilson coefficients 
C"f S to the known Wilson coefficients Cj MS correspond- 
ing to the operators Oi of basis I or II. 

To this end we first note that Eqs. © and jTOj hold 
for Rl-renormalized operators without contributions of 
evanescent operators, i.e., 



O 



RI 



T 10 Q'f 







(58) 



with 



T 



/1/5 1 
1/5 • 
3 
2 



3/10 
X3/10 -1 



(59) 
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where the matrix T can be read off from Eqs. ([9]) 
and (fT0|) . This equation holds for the operators Oi of 
basis I and II. Note that Q^-Qs are equal to Q' 5 -Q' s , see 
Eqs. (fTU|) . and therefore the respective sub-block in T is 
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proportional to the unit matrix. In the MS scheme, how- 
ever, finite contributions of evanescent operators modify 
Eq. ([SHI) to 

C>? s -T ij Q i f S = *t i G 1 (60) 
at the one-loop level. The coefficients Aij are given by 

At 4 = -— , Ati =0 for i^A (61) 

47T 

for operators Oi of basis I and 



la, . a, 
At 2 = 1 , Ai 4 = 

AU=0 for * 7^ 2,4 



(62) 



for operators Oi of basis II. 

Since the Hamiltonian is independent of the choice of 
operator basis, we have 



7 10 

■/MS/ x /^,/MS/ 



E CiM (qTm) - E ^(m) (of S M 



(63) 



Therefore using Eqs. (1501) and (Gi) = (Q p ) we can de- 
termine AT in 

cf%) = CpOO(T tf +AT«). (64) 
At one-loop order one finds [11] 

/ 



AT MS 



4?r 



(65) 



for the Wilson coefficients Cj corresponding to basis I 
and 



AT* 18 = 



47T 



for the Wilson coefficients C* 1 corresponding to basis II. 

In App. [B] we give an alternative way to determine 
AT using our results for an RI scheme derived in the 
next section. 



IV. CALCULATION AND RESULTS 

In the following we give the main results of this work. 
We first express the finite part of the MS[NDR] ampli- 
tudes defined in Sees. IIII CI and MI Dl in a compact and 
instructive way in Sec. lIV Al We then discuss the general 
matching procedure in Sec. IIVBI before giving the con- 
version matrices for d iffer ent R I/MOM and RI/SMOM 
schemes in Sees. II V Cl and ITvTJl 



A. A compact expression for the amplitudes 

The RI schemes are defined by applying projectors P n j 
in spinor, color, and flavor space to the renormalized am- 
putated Green's functions T^Of 1 ), see Eqs. ([55 ]) , ([47 )1 . 

and (51]), and thus to E^Of 8 ), see Eqs. (j3"5 |) and (gg j) . 



In this section we calculate r^ s (Of s ), r^ s (Of s ), and 
r|^(Op) as defined in Sees. iLTTCl and HITTfl The di- 
agrams are generated with the program QGRAF [53| . 
and the symbolic manipulations are carried out using 
FORM [5J|. We give results for the respective excep- 
tional momentum configuration of the RI/MOM scheme 
as well as the non-exceptional momentum configuration 
of the RI/SMOM schemes. We consider the operators 
Of 1 of basis I, II, and the chiral basis. 

We simplify the MS-renormalized amplitudes using the 
decomposition [3(| 

ItiHulp = g^ulp + gvpln - 9p, P lu + is a ^u P l a lb (67) 
and the Fierz identities 

W - 1*)\<iW - = W - 7fc)]«[*(l ~ 75)W 

-(g 2 /2)[7p(l-7s)M7 p (l-75)k, (68) 
[7p(1-7s)M7 p (1-75)]« = 

~[lp(l-l 5 )UY{l-ls)] k3 - (69) 

The resulting one-loop expressions are written as 

ifS(op) = rr c (o.) + g(T7rr e (0i) 
+ E ^T e (xf) 

k=q,Pi,P2 

+ TLlf e (^-) + Tf i rT"(Gij) , (70) 



( 66 ) ™ — 



pMS (0 M S) = ^l X G lr tree (Gi , 



47T 



(71) 
(72) 



where r^°°(Oj) denotes the insertion of the operator Oj 
of basis I, II or the chiral basis at tree level and the sum 
over repeated indices is implied. The operators Xj are 
obtained from the operators Oj by replacing the gamma 
structure 7 M ® 7 M with ft (g) ft/k 2 , i.e., 



X l = ( s a0 ~ 75)wb)("fc^(l - J 5 )d a )/q 2 



(73) 
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for the case of Qi and k — q and analog for all other 
Oj and k £ {q,Pi,p 2 }- Similarly the operators Yj are 
obtained from the operators Oj by replacing the gamma 
structure 7^ ® 7^ with ie^vapl^ ® l u PiP 2 lq 2 , i- e -' 

ii = ie MI y Q /3(s a 7 M (l - 75)wb) 

x(%f(l-75R)^/? 2 (74) 

for the case of O/i and analog for all other Oj. In 
Tabs. HHIIVI we give the coefficients T z with z £ 
{7, <f_,f v f 2 , e, Gi} at one loop for the chiral basis and 
the exceptional and non-exceptional momentum configu- 
ration. The corresponding results for basis I and II can 
be obtained from these tables using Eq. (|60p . The con- 
stant Go is defined as 

Go = (2/3)#'(l/3) - (2tt/3) 2 w 2.34391 , (75) 

where 'J (a:) is the digamma function [551 ] . 

In order to determine the conversion factors from the 
RI schemes to the MSfNDR] scheme we need to study 
projected Green's functions 

A nst = P nt T^ (Op) , (76) 

where n G {2,4,4p} and both P nt and T^ s have open 
spinor, color, and flavor indices. The action of P„t on 
rj^ s for n £ {4, 4p} is thus given by 

A _ pa/3t8;ijkl;f r MS //lMS) / 77 \ 

JVnst — -r„i 1 n;a/3j6;ijkl;f\ u s ) \" J 

with spinor indices a, /3,7, <5, color indices i,j,k,l, and 
flavor index /, as shown in Fig. [31 Similarly 

where [i is a Lorentz index and a £ {1, . . . , A^? — 1} enu- 
merates the generators of the SU(iV c ) algebra, see Fig. [5] 
The summation over repeated indices is implied. The re- 
sulting expression A„ st naturally depends on the choice 
of the projectors P n t- In the following we make certain 
assumptions about the structure of the projectors, i.e., 
we discuss different classes of projectors. For each class 
of projectors we can then give a very compact form of 
the amputated Green's function that can, without loss of 
generality, be used to calculate all projections A of the 
respective class. 

First we restrict the discussion to projectors which con- 
tain no external momentum, except for at most the ex- 
ternal momentum 

Q = Pi - P2 ■ 

We also do not consider projectors that contain vectors 
such as 6 pi, 6^2, <5p3, or S^a which break Lorentz sym- 
metry explicitly. The corresponding class of projectors 
shall be denoted by Pi, We will also study the sub-class 
P7m £ pi f projectors that additionally do not contain 
the external momentum q. 



Let us focus on the amputated Green's functions with 
four external quarks and examine the projection of a 
spinor structure such as fliiGSifia under projectors of class 
Pi, i.e., 

where the two factors belong to the two fermion lines, 
and the combination fli ® contains strings of gamma 
matrices without remaining open Lorentz indices. Then 
the ansatz 

P% M ni ® 7„fi 2 ] = hg^ + hq^/q 2 (79) 

with Lorentz indices fi and v is justified by the Lorentz- 
transformation properties of both sides since Pi does not 
contain any Lorentz vectors apart from the momentum 
q. The coefficients t\ and t 2 can be determined by con- 
traction with g^v and q^qv, respectively. We find 

P% M fii<»7^2] = 

+ |^ (-PtbpSll ® 1^2} + AP*[$li ® ^ 2 }/q 2 ^j 

+ 9j f [Pthpttl ® 7^2] - P'[0fii ® ^ 2 }/q 2 ^j . (80) 

Therefore not all of the spinor structures that appear in 
the general MS amplitude of Eq. ([70)1 are independent 
under projection with Pi, i.e., 

p?[rr°(Af )] = pt[TT c {xf)\ = \ptmOj)] , 

Pi[T^(Y j )}=0 (81) 

with n £ {4,4p}. The resulting projected amputated 
Green's functions for n = 4 up to one-loop order can be 
written as 

p^rY s (oY s )] = p*[r t ™{o i )} 

+ Tjp^[r* rec (^)] + Tf 'P'tr^Gi)]) . (82) 

If we further restrict the discussion to projectors of the 
sub-class P 7 " C Pi, we find 

P^ [TT c (Gi)} = \p~<» [TT C {Q P )\ (83) 

with n £ {4, Ap} using the same arguments as in Eq. (|79[) 
without the term proportional to t 2 . Therefore, for pro- 
jectors of class P 7m and n = 4we can write 

PMrf^op)] = P ^ rrm] + ^((r 7 + It J 

+ i T & + i T « + -^f^P^mOj)]) (84) 

with r4 rco (Qp) = TjT^ ee (Oj), where we can read off the 
coefficients Tj from Eq. P^)l for basis I and II or from 
Eq. (|2"9")l for the chiral basis. 
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1.1) C(~ ±] t i2i +41og(2)+iV c + 1 i--l)-H^21+121og(2) +1 I--7 (-2iV c -2)£ 

2.2) e(-l^+AT c+ i)-i«l + £ _2iV c C 

2.3) (41og(2)-l)£+121og(2)-7 -2£ 

3.2) (41og(2)-l)£+121og(2)-7 -2£ 

3.3) Z(-^+N c +±)-™%fl+£ -2Nc£ 

c c\ tl 4 1og(2) | ^ 5 \ 41og(2) 5 t / 81og(2) 9 jy , 8 \ 81og(2) 4 n n 

°>°J ?l W c hJVc 3A^,< 3JVc 3iV^ ?l 3A^ ^ Vc +3A^J ~3W C 3JV^ U U 

5,6) (Il2l(21 + i^ + ii2i<21 + i (»iaei_|) €+ »iaei + | 

6,5) («2fW_I)£ + ±i2|M_| (8]M2)_|) C+ 8iM2) + | 

a p\ tt 41og(2) , o » r _ 5 \ 41og(2) , » r 5 f/ 8 81og(2) ^ 81og(2) 4 

U ' U ^ M, 3JV,. tZiVc gjVeJ 3JV. "t^JVc 3jVf ^VANr 3N r J 3N^ 3N r U U 



TABLE II. One- loop coefficients T 7 , T' 4 !, Y*, and T e for the exceptional momentum configuration and the chiral basis. The 
coefficient Y£? = for the exceptional momentum configuration. The coefficients for Q7 and Qs are identical to the coefficients 
for Q$ and Qa. The remaining matrix elements not given here are zero. 

V2_±l i.i i.i ' a ZJi Zli 

(1.1) e (_2£|^ + ^_2|0„41ogi2) +41og(2) „2CaA t _2Co + (_2| £ _|^ ^^ID 5 

-L9N 1 ^| , CojVe , Co 12 log(2) + Mo + 4 ' _2) . M<L + 4 

+ iiY = _ l^J H 3 h ~ WT~ 3 3 ■ i ' 3 3 

+121og(2)-2JV c +jf-9 

(2.2) ? (-2C ^ + ^_iiM2)+27V c - 1 I-) _20^_2iM + 4^ 1^)+^ 

, C iV c 12 1og(2) n * r r 9 
+ 3 JV C ZJVc_I "JVc 

(2.3) e(41og(2)-2|o)+^+121og(2)-9 -2|o-f + § (2Sl_|) €+ | 

(3.2) ? (41og(2)-2fo)+^+121og(2)-9 -H|o_H+| (3Qi_|) f+ g Q 

(3.3) e(-^^ + ^- 1 ^ 21 +2iV c - 1 i-) _2CoA t „2iV ££+ 4| £ e (2£ ^„2A t ) + 4| £ Q 

, CgNc 12 1og(2) 9^- , 9 
~r 3 JV C ZJVc+ JV c 

^ ^ 2C iV c , 7C 41og(2) v §_\ 2C iV c , 2Cp , t l 41og(2) 2iV c t( 2CgN c 2Cp 2iV c Cpg 2C iVe 

°y SV 3 6iV c 3JV C " rz,1,c 3AT C J 3 3JV C ' ? I 3N C 3 V 3 3iV c 3 3JV C 3 

, C iV c , 7Cp 41og(2) 0N , 1 l \ 41og(2) 4iV c 2 , 2_ \ , 4iV. _ _J_ , Cfc 

3 6iV c 3iV c ZJVc+ 3iVl + JV^J 3JV C 3~~ JV^ 3iV c / 3 3AT C ^ iV c 

(5,6) ? (_^+ii^ + |)- 3 fi + ii£|(2) + | (iMM-D^+fi^+l o ai-Si 

( 6 ,5) ? (_2| a + ii2l(H) + |) + C|l + i]£l(2)_W _3C iH .(l]aei_l) f+ l]aei + 2 (2|0_| )c+ 4 _2Cp 

Cfi fil CpiVc I 7Cp 41og(2) . r, y 8\ 2C , f / 1 41og(2) \ 41og(2) */_2 2Co\ 4_ t( C N c Cg \ 

q,y 2 t 6JVc 3JVc T^J'c 3iVc ; 3JV C """^ \ JV C 3N C J 3N C >> \ 3JV C 3iV c ^ 3N C >\ 3 3N C J 

3C N C i 7C 41og(2) , H Ar | 1 2_ gQ^c i Cp 

2 ^SiVc 3iVc ~ l ~^ JVc+ 3iV e Na 3 J_ N a 

TABLE III. One-loop coefficients T 7 , T^i, T^2, T^, and T £ for the non-exceptional momentum configuration and the chiral 
basis. The coefficient T^j = for the exceptional momentum configuration. The coefficients for Q-j and Qs are identical to 
the coefficients for Qs and Qe- The remaining matrix elements not given here are zero. 



*Y*"l "v"* G x G i ry-^ G i r y r G\ r *fG± G i 
1 1 x 2 x 3 X 5 1 6 1 7 

-f -§ 

TABLE IV. One-loop coefficients T Gl for the chiral ba- 
sis. The results are identical for the exceptional and non- 
exceptional momentum configuration. 



B. The matching procedure 

In the following sections we give the one-loop con- 
version coefficients from different RI schemes to the 
MS[NDR] scheme. We first explain the general procedure 
of calculating the conversion coefficients and then give re- 
sults for explicit RI/MOM and RI/SMOM schemes. 

At one loop the conversion from RI schemes to the 



MS[NDR] scheme in terms of renormalized amputated 
Green's functions reads 

OT^(op) = r^of 1 ) + Aaf^ri(o? ) 

+ Ac RI^MS r , (G RI ) (85) 

with wave function conversion factor C v — (CV) 2 for 

n G {4,4p} and C v = C y q {Z^ / Z y A ) 1 / 2 for n = 2, see 
Eq. (|53|) . In order to define the RI scheme, and hence to 
determine the coefficients Ao^ hMS and Ac? 1 "^ 8 , we 
apply projectors P n i to Eq. (|85[) and use the RI condi- 
tions of Eqs. (|38]) and gT]) for the RI/SMOM 2 schemes 
and of Eqs. (J35J) and (J3TJ for the RI/SMOM schemes, see 
SecHHDj 

From Eqs. dTTJ) and $TZ§ and the RI conditions of 
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Eqs. (1471) and (|51j) we can already conclude that 



C. The RI/MOM and RI/SMOM(7 (1 ,y) schemes 



_ " 47T 



(86) 



is the only possible result at one loop. For convenience 
we provide a specific projector to determine the mixing 
of the two-quark operator G\. In an RI/SMOM 2 scheme, 
we can use 



(87) 



where A° is a generator of the SU(iV c ) group algebra and 
the convention of open indices is given in Fig. [5] In an 
RI/SMOM scheme, we can use 



Pav,G x = 5il8k 3 {"f P ')pu{ln)5 1 ■ 

At higher loops the additional operators [3£ 



(88) 



We first give the conversion coefficients for RI schemes 
that either (i) use only projectors of class P 1 * or (ii) use 
the exceptional momentum configuration and projectors 
of class Pt. The coefficients Ac RI ^ MS are already deter- 
mined in Eq. (|86[) . and it remains to obtain the coeffi- 



cients Ai 



,ri-s.ms 



by applying projectors Pik to the am- 



MS 



putated Green's function T\ 

In case (i) the amputated Green's function sim 
plifics to 

PM(c«) 2 rf3(oMS)] 



P^[T% ee (Oi)] 



4 y 



4 y 



It* 

4 y 



3„c 

4 1 



x pT"[r^ roo (Oj)] 



(93) 



G 2 = ■r T s{ J D M £»' i , J D ,/ }7 I/ (l-75)d, 

G 3 = -^aDpDuDxS^Q. - 7&)d, 
«ff 



(89) 



with S^"* = j^"f^j x — ry x j lJ -fV and possibly other non- 
gauge-invariant two-quark operators N m mix with the 
(8,1) operators. We therefore need to provide additional 
renormalization conditions at higher loops in order to 
define the RI schemes for the (8, 1) operators uniquely, 
i.e., we need to specify additional projectors. The spe- 
cific choice of projectors to determine the mixing with 
the operator G2 and additional operators which occur at 
higher loops is not relevant for the one-loop conversion 
factors given in this work. At one loop the conversion 
factors can thus be given without specifying the details 
of higher-loop contributions. 

In the following we use the RI'/MOM, RI/SMOM, and 
RI/SMOM wave function renormalizations defined in 

Refs. O Ea O- 

The corresponding conversion factors 
up to one-loop order 



1 + —A y 

4tt q 



are given by 



a RI'/MOM 

9 



^RI/SMOM 



2\ N n 



(90) 



(91) 



for the RI'/MOM and RI/SMOM wave function renor- 
malization schemes as well as 



RI/SMOM^ _ if 1 

- 2 r c Nr 



for the RI/SMOM scheme. 



-i+i(3-c„: 



(92) 



under projectors P^k £ P 7m , see Eq. 
be compared to Eq. ([85]). i.e., 



This should 



pv»{(cy) 2 r™ s (o™ s )} = p^irKof 1 )] 



( 



Aa 



RI^MS 



§Ac™r>^[rl(Of)], (94) 



where we used Eqs. (|83j) . If we impose the RI condition 
of Eq. (J3HJ) and insert Eq. USB}, we find 



Aa 



RI->MS 



4tt V 



T7 



I(T^i 
4 



ft 2 1 t^ 



2A y q l 



(95) 



with identity matrix 1. 

In case (ii) we have 
simplifies to 



0, and therefore Eq. 



ps*[rf s (of s )] = pi*[r* rcc (o 2 )] + — f(T 

47T V 



4 v 



-Tfnp^[r* rcc (o 3 )] + 1 j^[r* roe (Gi) 



(96) 



One can read off the conversion matrix Aa RI_ *' by com- 
paring Eq. dSSJ to Eq. ([96j) with 



^ RI->MS _2±('T1 

4tt V 



It^i + It^ + 2A v n 1 
4 4 9 



(97) 



In both cases the conversion coefficients are unique up 
to the definition of the RI quark wave function renormal- 
ization A^. In other words, in these cases the details of 
the projectors do not matter, as long as one uses a suffi- 
cient number of independent projectors to determine all 
elements of the conversion matrices or Z-factors, respec- 
tively. 

The choice of the RI wave function renormalization 
scheme only affects the diagonal elements of Aa RI ~ ! ' MS . 
One can also combine different RI wave function renor- 
malization schemes with different RI schemes for the op- 
erators Oi in a straightforward way using Eq. (|95l) . 
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We call the RI scheme of case (i) RI/SM0M(7 M , y) 
scheme, where y — d corresponds to the RI/SMOM 
wave function renormalization scheme and y = 7 M corre- 
sponds to the RI/SMOM^ wave function renormaliza- 
tion scheme. The name RI/SMOM(7 jU , y) reflects that 
we restrict the choice of projectors to the class P 7 ". 

The scheme corresponding to case (ii) is the RI/MOM 
scheme, and we give results only for the RI'/MOM wave 
function renormalization. 

In Tabs. lVHVlIl we give the values for Ar RI ^ MS , which 



is defined in Eq. (j56j) . for the RI/MOM scheme with 
RI'/MOM wave function renormalization and operators 
of basis I, II, and the chiral basis. In Tabs. lVnll and llXl wc 
present the results for Ar RI ^ MS for the RI/SMOM(7 M , jjf) 
and the RI/SM0M(7 M ,7 Al ) scheme and operators of the 
chiral basis. 

We would like to point out that the result of Tab. [V] 
for iV c = 3 and £ = or £ = 1 agrees with the result 
of Ref. [32j for the conversion from the RI/MOM scheme 
to the MS[NDR] scheme. In Ref. [30[ the current-current 
contributions in the RI/MOM scheme have been consid- 
ered. Our result for these contributions agrees with the 
one of Ref. [30( , which can be seen by using the results of 
Eq. (|60[) and Tab. Inland combining them with the quark 
wave function renormalization constant of Eq. (5.3) of 
Ref. [11. 

Furthermore, the result for the (27, 1) operator Q \ in 
Tabs. [Vm] and E] agrees with the result of Ref. [28| 
for the conversion from the RI/SMOM(7 AI , and the 

RI/SMOM(7 M ,7 M ) scheme to the MS[NDR] scheme for 
the (W+AAHAS = 2) operator. 



D. Results for the RI/SMOM (j^, y) schemes 

In the following we discuss the non-exceptional mo- 
mentum configuration and RI/SMOM schemes defined 
by projectors of the class P$ . This allows for the defini- 
tion of new and independent RI schemes which we call 
RI/SMOM(;|, y) schemes, where y denotes the choice of 
the wave function renormalization scheme as in the pre- 
vious section. We give results only for the chiral ba- 
sis which does not contain linearly dependent operators. 
Since mixing only occurs within the blocks of the (27, 1), 
(8,1), and (8,8) operators, we renormalize each block 
separately. 

We first note that for a scheme with projectors of the 
class P$ the conversion matrices cannot be simply read 
off from Tabs. HTMlVI due to the nonzero contribution of 
T0. This means that the schemes defined in this sec- 
tion, unlike the schemes discussed in the previous sec- 
tion IIVCI depend on the specific choice of the projec- 
tors. However, any non-degenerate linear transformation 
of a given set of projectors leaves the conversion matrix 
of Eq. (|3"4"|) invariant. We give the projectors used to de- 
fine the RI/SMOM(rf, y) schemes in the following. The 
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(41og(2)-|)C+121og(2)-f 
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9N C 
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Hog (2) jy, 



(21og(2) + i)e+21og(2)+2 



21og(2) 2_ 

N c N c 



(21og 

(2 log 
(2 log 



(2)-l)e+21og(2)+3^-2 

2 1og(2) \ 21og(2) , » r _J2 5 

N c I N c ~r* I, c Nc 3 
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N a 
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(4 log 
(4 log 
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12 1og(2) 7 
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(2)-§)£+121og(2)-7 
(2)-§)£+121og(2)-7 



og(2) N c , 2 \ 12 1og(2) . 7 
N r 2 ' N r . J Nc ' Nc 



-0.43926 
1.31777 
1.31777 

-0.43926 
0.07407 

-0.22222 
0.07407 

-0.22222 

-0.29111 
0.87332 
0.14815 

-0.44444 
1.87332 

-2.10592 

0.55556 
-1.66667 
-1.12876 

3.38629 

0.55556 
-1.66667 
-0.05815 

9.20457 
-1.12876 

3.38629 
-0.61371 
10.87124 
-0.07407 

0.22222 
-0.07407 

0.22222 
-0.43926 

1.31777 

1.31777 
-0.43926 



TABLE V. One-loop conversion matrix Ar RI ^ MS for basis I 
in the RI/MOM scheme with RI'/MOM wave function renor- 
malization. We also give numerical values Ar* for N c = 3 
and f = 0. 



(i,3) 


l J I 47T 


Ar h/ii 


(2,3) 


5 

9iV e 


0.18519 


(2,4) 


5 
9 


-0.55556 


(2,5) 


5 
9iV e 


0.18519 


(2,6) 


5 
9 


-0.55556 



TABLE VI. One-loop conversion matrix Ar RI ^ MS for basis 
II in the RI/MOM scheme with RI'/MOM wave function 
renormalization. All elements not given here are identical 
to Tab. [V] 
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hJ) 

mT 

2,2) 
2,3) 
3,2) 
3,3) 

3,5) 
3,6) 
5,5) 
5,6) 
6,2) 
6,3) 
6,5) 
6,6) 
7, 7) 
7, 8) 
8,7) 



Ar 



RI->MS / a 



/a* 
I iir 



^(-^+4108(2)-^+^-!] 



-H±2£i2>+121og(2)+^-7 

41og(2) N c , 2 \ 12 1og(2) 
2 ' N c 



(41og(2)-f)C+121og(2)-7 



07 



41og(2)-|)C+121og(2)+3 W --^ 

41og(2) JV 2 \ 12 1og(2) 67 
JVc 2 ' N C J N c 9N C 



f-l 4 1og(2) JV C ^ , ^,,^ s ^ ; 

> V JVc 2 ' N r ' K - ' 

_2 
3 
2 

9 JV C 
_ 2 
9 



21og(2) JV C \ 21og(2) 
JV C 2 / JV C 



2 
JV,, 



(21og(2) + i)C+21og(2)+2 

_5 10 

JV C 3 

3JV C 

(21og(2)-l)e+21og(2)+4--2 



2 1og(2) 
JVc 

c/ 21og(2) JV, 
^ V JV C 



21og(2) 



JV, 
21og(2) 
JVc 



3JV C 

+4jV c 



2 

' JV, 



■2 

N 



(21og(2) + ±)M-21og(2)+2 
(21og(2)-l)£+21og(2)-2 



.8) £(N C - 



2 1og(2) N 
JVc / 



21og(2) 



JVc 



0.87851 

-0.43926 
1.31777 
1.09554 

-0.95777 



0.07407 
-0.22222 
-1.12876 

3.38629 
-1.66667 
-3.88889 
-0.05815 

9.20457 
-1.12876 

3.38629 
-0.61371 
10.87124 



TABLE VII. One-loop conversion matrix Ar R1 ^ Mb for the 
chiral basis in the RI/MOM scheme with RI'/MOM wave 
function renormalization. We also give numerical values Ar* 
for N c = 3 and £ = 0. 



projectors will be expressed in terms of 

VV±AA,f 

P {1)J , = dff'dijdki 

VV±AA,f 
P (2)J> = 6 ff' 6 il d kj 

VV±AA4 . - 

P {3)J , = SffOuSkj 

x [{0p-y(0Sa ± (^75)^7 {tl5)s a ]/q 2 , 

x [{i)M{i)s a ± (^W^sM/V , (98) 



P, 



VV±AA,, 
(4),/' 



where all indices are as shown in Fig.H 

The (27, 1) operator Q[ only mixes with itself, so we 
only need to provide a single projector P? 27 j>. We adopt 

the RI/SMOM^,0) and RI/SMOM(<£,7 M ) schemes de- 
fined in Ref. [28| so that for the (27, 1) operator we use 



P, 



P, 



VV+AA.f 
(!),» 

(27,1), 1 6 4iV c (iV c + 1) 



(99) 



The (8, 8) operators Q' 7 and Q' s also only mix with them- 
selves under renormalization, and therefore we have to 
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TABLE VIII. One-loop conversion matrix Ar RI ^ MS for the 
chiral basis in the RI/SMOM(7 M , scheme. We also give 
numerical values Ar* for N c — 3 and £ = 0. 



provide two projectors. We choose 
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The one-loop mixing coefficients for the (8, 1) operators 
shall be determined using the projectors 
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TABLE IX. One- loop conversion matrix Ar 



for the chi- 



ral basis in the RI/SMOM(7 M , 7 M ) scheme. We also give nu- 
merical values Ar* for N c = 3 and <; = 0. 
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C(-%+21og(2) + i)-2Qi+21og(2) 

>V 2 o \ / 1 2 / 2 o v / 

+2 


-0.12957 


(6,2) 


5 10 

JV C 3 


-1.66667 


(6,3) 


IS Cj 

3JV C 


-3.88889 


(6,5) 


(21og(2)-I)C+21og(2)+3^-2 


-0.05815 


(6,6) 


3C JV C , 3C c(N c 2 1og(2)\ 
2 1 2JV C H 2 JV C / 

21 ; s(2) +5jv c ^ § 


2.49561 


(7 7) 


f /C ^c 2 1og(2) JV C \ . 3C 2 1og(2) 
?V 2 JV C 2 J 1 2JV C JV C 


2 70986 


(7, 8) 


£(-%+2 log(2)+|) -2§a+2 log(2) 

^ V 2 o \ / 1 2 / 2 o v / 

+2 


-0.12957 


(8,7) 


(21og(2)-i)C+21og(2)-2 


-0.61371 


(8,8) 


3C JV C | 3C | cz-JVc 2 1og(2)\ 
2 1 2JV C M 2 JV C / 
21og(2) , K Ar 3 
JV r + 0JVc JVc 


4.16227 



TABLE X. One-loop conversion matrix Ar RI ^ MS for the chi- 
ral basis in the RI/SMOM(^, 7 M ) scheme. We also give nu- 
merical values Ar* for N c = 3 and £ = 0. 



where 



j = 2,3,5,6. 



(103) 



Other choices of projectors within class Pt, which are 
not linear combinations of the projectors given above, 
are possible and might lead to smaller conversion factors 
and a better convergence of the perturbative expansion. 
The reader can easily obtain the respective projections 
using the results which we provide in Tabs. UTrllVI 

In Tabs. |X] and IXII we provide the resulting one-loop 
conversion matrices Ar RI ^ MS for the RI/SMOM(^, 7 M ) 
and RI/SMOM(^) schemes. We note that the result 
for the (27, 1) operator agrees with the result of Ref. [28[ 
for the conversion of the (VV+AA)-(A5 = 2) opera- 
tor. In Tabs. IXTlrlXIIII we give the individual conversion 
matrices Aa RI ^ MS for the (8, 1) operators in the same 
schemes. 



E. Discussion 

In the previous sections we provided results for the con- 
version from different RI/MOM and RI/SMOM schemes 
to the MS[NDR] scheme. In general the different 
schemes have a different magnitude of one-loop correc- 
tions Ar RI_> and a different rate of convergence of the 
perturbative expansion of Ar RI_>,MS in a s . In the context 
of lattice QCD and non-perturbative renormalization it 
is thus useful to have several RI schemes to choose from 
in order to estimate the effects of missing higher-order 
terms in the perturbative expansion. In this work we 
give four different RI/SMOM schemes for each operator 
of the chiral AS = 1 operator basis. 

Since we only give one-loop results, we cannot estimate 
the rate of convergence of the different schemes presented 
in this paper. In some cases where higher-order results 
in RI/SMOM schemes are known, such as the conver- 
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h J) 



Ar 



RI->MS / a 



I iir 



1,1) 

2,2) 

2,3) 
3,2) 
3,3) 

3,5) 
3,6) 
5,5) 

5,6) 

6,2) 
6,3) 
6,5) 
6,6) 

7, 7) 

7, 8) 

8,7) 
8,8) 



12 1og(2) 
N r 



e(f-Co-«+41og(2)> 
+121og(2)+^-9 

t ( 4C N c . Cp 6C 41og(2) 4iV c , 
H 5 T « c 5 w. 

12 1og(2) 8iV c 9_ 

N c ' 5 ' N a 
t( 6C N c 9Cq 



Nc 
12 
5 



+41og(2)-S^+| 



+121og(2)+if^-f 



6Af c 



e( _6C^ + 4Co +41og(2) + ^ 
12Af c 



+ 121og(2)- 



t ( 9C N C . Cp , 6C 
S V 5 T » C T 5 

18iV c 
5 



41og(2) , 9JV, 



12 1og(2) 



JV C 



i 85 , 26 
~t~ Q AT - I 1 K 



2 1og(2) 

iV e 2JV, 



3C-H 
2JV, 



2 1og(2) 



-!)- 

2 

9iV c 
_ 2 
9 

_C„ 

2N C 
2_ 

JVc 

e(-^+21og(2)+i)-«fa+21og(2) 
+2 

_5 10 

iV c 3 

^2--5 

3JV C J 

(21og(2)-i)C+21og(2)+3^-2 

CgN a , Cp 21og(2) » r l\ 

2 ~ r 2JV c AT C Tiv c 2iV c 7 

3C iV c , 3C 21og(2) . y _^ 5 

2 ~ r 2W c iV c Tiiv c ^ 3 

t( Cp 2 1og(2) 1 \ , 3C 2 1og(2) 

5 V 2JV C JV C 2JV C / ~ 2N a JV C 
2 



e(-^+21og(2) + i)-^+21og(2) 
+2 



(21og(2)-i)C+21og(2)-2 
f/ C n JVe i gp 21og(2) | Ar 



2 

3C iV c 



2iV e ]V C 
3C 21og(2) 



2JY, 



iVc 



+JV C 
+4iV c 



i 

" 2JV e 
2 

' iVc 



-0.45482 
2.62741 

4.91777 
-4.50446 
-8.69111 

0.07407 
-0.22222 
0.04319 

-0.12957 

-1.66667 
-3.88889 
-0.05815 
-0.17106 

0.04319 

-0.12957 

-0.61371 
1.49561 



TABLE XI. One-loop conversion matrix Ar RI_>MS for the chi- 
ral basis in the RI/SMOM(^, scheme. We also give numer- 
ical values Ar* for N c = 3 and £ = 0. 



A RI-»MS / a. 



(2,2) £(±^ + 



Cp 

2.V, 



(iCp 
5 



12 1og(2) . 13iV c . 8 

N c ' 5 ' Nc ' 



41og(2) 13iV c 
iV c 10 

12 



1 

" 2JV, 



(2,3) e(^-^+41og(2)-^ + |)+121og(2)+lfi 



(3,2) e(-^+ 4 § a +41og(2)+M £ -f)+121og(2)- 



j-j.v. 

5 



(o o\ <;/ 13C JV e , Cp , 6Cp 41og(2) 13iV e , 1 en 

\°1 °) <i\ 10 T 2Nc ^5 iV c 10 2iV c 5 / 

121og(2) 13iV c j 8 ,12 

N 5 ' N c ' 5 

I'c r\ <-( C Q N c 21og(2) JV C \ , 3C Q 
1°>°J ?( 2 JV7 2~J + 2JV C 

(5,6) 5(-^+21og(2) + i)-5fi+21og(2)+2 
(6,5) (21og(2)-±)£+21og(2)-2 



2 1og(2) 



Nc + N -~4- c 



(6,6) 



3C iV c | 3C, 



+ 



2 -V, 



21og(2) ' 
JV C , 



2 1og(2) 



-5jV c 



3 



TABLE XII. One-loop conversion matrix Aa RI ^ MS for the 
(8,1) operators of the chiral basis in the RI/SMOM(^, 7 M ) 
scheme. 



I ■ -\ a RI-KM3 / a. 



, 9 12 

< /V. K 



C9 9^ cf iCgNc , Cp 6C 41og(2) 4JV C , 6\ 
l/l Z ) ?( 5 h JV7~~ 2V C 5 h sJ JV C 

+ ■ " 

(2,3) e( 
(3,2) e( 

(3,3) 



121og(2) i 8JV, 



+- 



12JV C 



JV C 5 

6CoiV e 9C i ^ , ? ^o^ 6JV C 

5 5 

53 
5 

6C iVe | 4C j 41 g(2) | 6N " 9 > !■> I. .,./•■>! i-'.v. 

27 



+41og(2)-^ + f) + 121og(2) + f 



+121og(2)-ifk 



mCgj 6Cq 41og(2) 9JV e 
" "r" JV C • 5 JV C "r" 5 



6> 



121og(2) 
JV C 



18JV e i 9 I 12 
5 ' JV C "•" 5 



.2JVC 



2 1og(2) J_ 

JV„ 2JV C 



3C 2 1og(2) 
' 2N a JV C 



(5.5) e( 

(5.6) e(-^+21og(2) + I)-3fo+21og(2)+2 
(6,5) (21og(2)-i)e+21og(2)-2 

/ fi fi-l d-f Cp JV C , Cp 21og(2) . N 1 \ 3C JV e , 3C 

\v,v) c;y 2 -|- 2JVc Nc -t-Jv c 2JV(= j -f 



2 ' 2JV e JV C 



2JV C 



JYc_ 



JYc_ 



TABLE XIII. One-loop conversion matrix Ao B 



for the 



(8,1) operators of the chiral basis in the RI/SMOM(^, i 
scheme. 



sion relation for light quark mass determinations, the 
convergence in the RI/SMOM schemes is significantly 
faster than the convergence in the traditional RI/MOM 
schemes. 

It is, however, interesting to compare the magnitude 
of the one-loop corrections Ar RI ^ MS of different RI 
schemes. In Tab. IXIVI we give the spectral matrix norm 
and the maximum norm for the conversion matrices Ar* 
at N c — 3 and £ = for the different blocks of operators 
in the chiral basis. The spectral norm of matrix M is de- 
fined as the maximum singular value of the matrix M and 
the maximum norm of matrix M is defined as the max- 
imum absolute value of the matrix elements. Therefore 
the spectral norm gives a good estimate of the general 
magnitude of the one-loop corrections since it is invari- 
ant under a change of operator basis. The maximum 
norm gives a good estimate of especially large individual 
mixing coefficients. 

One observes that for the (27, 1) operator the one- 
loop coefficients for the RI/MOM scheme with RI'/MOM 
wave function renormalization are of order 1, while the 
RI/SMOM results vary from order 1/10 to order 1. The 
RI/SMOM(a;, y) schemes with x — y have especially 
small one-loop coefficients. In this case the spectral 
norm is, of course, equal to the maximum norm. For 
the (8,1) operators the spectral norm and the maxi- 
mum norm for the RI/MOM, RI'/MOM scheme and the 
RI/SMOM(^,y) schemes is approximately twice the size 
of the RI/SMOM(7 M , y) schemes. In the case of the 
(8, 8) operators the maximum norm as well as the spec- 
tral norm in the RI/MOM, RI'/MOM scheme is signifi- 
cantly larger than the respective norm in the RI/SMOM 
schemes. The difference is especially pronounced for the 
RI/SMOM((j, scheme. 

We conclude that the magnitude of the one-loop con- 
tributions to Ar RI ^ MS varies significantly between the 
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Scheme 


(L,R) 


l|Ar*/£l|. 


IIArVSflk 


RJ/MOM, RI'/MOM (27,1) 


0.87851 




RI/SMOM( 7m ,^) 


(27,1) 


2.45482 




RI/SMOM( 7m , 7m ) 


(27,1) 


0.21184 




RI/SMOM(^, fj 


(27,1) 


0.45482 




RI/SMOM(^, 7m ) 


(27,1) 


2.21184 




RI/MOM, RI'/MOM 


(8,1) 


10.61021 


9.20457 


RI/SMOM( 7m ,^) 


(8,1) 


5.34301 


3.88889 


RI/SMOM( 7m , 7m ) 


(8,1) 


5.21189 


3.88889 


RI/SMOM(^, fj 


(8,1) 


12.03242 


8.69111 


RI/SMOM(^, 7m ) 


(8,1) 


11.17384 


6.02444 


RI/MOM, RI'/MOM 


(8,8) 


11.42390 


10.87124 


RI/SMOM( 7m ,^) 


(8,8) 


3.23249 


2.62348 


RI/SMOM( 7m , 7m ) 


(8,8) 


4.77841 


4.49561 


RI/SMOM(^, <f) 


(8,8) 


1.62236 


1.49561 


RI/SMOM(^, 7m ) 


(8,8) 


4.26036 


4.16227 



TABLE XIV. Magnitude of the one-loop conversion matri- 
ces Ar* / for operators in the chiral basis and N c — 3, 
£ = 0. We give the spectral norm ||Ar*/^-|| 3 as well as 
the maximum norm 1 1 Ar*/|^-| |oo for operators in the (L,R) 
representation of SU(3)z, <g> SU(3)b. 



different schemes, and we expect that the same holds for 
the rate of convergence of the perturbative expansion of 
Ar ra^MS in Us Note that Tab. IXTVl only gives the nu- 
merical values for Landau gauge and that the qualitative 
analysis of this section may be different for other gauges. 



V. SUMMARY AND CONCLUSION 

Physical processes that change the strangeness by one 
unit play an important role in the field of flavor phe- 
nomenology. These processes can be studied in lattice 
simulations using an effective AS* = 1 Hamiltonian of 
electroweak interactions that is formulated in terms of 
AS = 1 flavor-changing four-quark operators. In order 
to renormalize these operators non-perturbatively and to 
convert measured matrix elements to the MS scheme it 
is necessary to define renormalization schemes that are 
independent of the specific regulator. In this work we 
define RI/SMOM renormalization schemes for different 
AS = 1 operator bases and provide one-loop matching 
factors to the MS[NDR] scheme. 

Since the different RI schemes project out different 
components of the perturbative series, the variety of RI 
schemes can be used to estimate the effects of higher- 
order corrections to the conversion matrices which are 
presently unknown. In this work we define four different 
RI/SMOM schemes for each AS — 1 operator. We also 
provide a compact expression for the finite one-loop am- 
plitudes that can be used by the reader to define further 
RI/SMOM schemes in a straightforward way. 

The numerical size of the one-loop contributions to 



the matching factors is discussed briefly for different RI 
schemes in the Landau gauge. We find that their magni- 
tude varies significantly between the different RI/MOM 
and RI/SMOM schemes. 

In future work two-loop conversion factors will be cal- 
culated which will further reduce the systematic un- 
certainties in lattice calculations involving the effective 
AS* = 1 Hamiltonian of electroweak interactions. 
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Appendix A: Flavor and isospin decomposition of 
four-quark operators 

In this appendix we give the flavor and isospin decom- 
position of the four-quark operators defined in Eqs. @- 
©. We proceed along the lines of App. B of Ref. [2j, 
carefully avoiding the use of Fierz transformations. We 
can then convert operators Qi to operators Qi including 
the evanescent operators of Eq. © . 

1. Left-left operators 

The left- left operators Q\, . . . , Q4, Qg, Qiq can be writ- 
ten as 

Qll = (TllTM ® <&M ® Ql) (Al) 

with left-handed quark fields q\ and flavor indices i, j, 
k, I. We identify q\ with an up quark, q\ with a down 
quark, and q\ with a strange quark. The color and 
spinor contractions are denoted by ®. The individual 
quark fields transform in the fundamental representation 
of SU(3) L , i.e., under V G SU(3) L we find 

Q'll = (Vi) io (Vl) jb (T LL )%V ek V < a(q i L ® q k L )(qi ® q l L ) 
= {T' LL )l{q i L ®qi){<? L ®q l L ) (A2) 

with 

(TL)t = (V^ a (V^ jb (TLL) a c b d V ck V dl . (A3) 

This transformation corresponds to the 81-dimensional 
representation (3 (g> 3) ® (3 (g> 3) of SU(3)/,. We can de- 
compose the tensor of this 81-dimensional representation 
as 

(T LL )% = (T LL m + (T LL )l$ 
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(A4) 



+ (T LL )f k % + (T LL )^, 

where for a general tensor we have 2tuj\ = ty + tji 
and 2t[j 3 ] = tij — tji. It is straightforward to check that 
the respective subspaces are invariant under Eq. (| A3|) and 
that their dimensionality is given by 



(lb) 



(All) 



dim 
dim 



(rp 



= 36 . 
= 18. 



dim 
dim 



(rp \ 
(rp 



= 18. 
= 9. 



(A5) 
(A6) 



Since Qll is symmetric under simultaneous exchange of 
(i, k) -H> (j, I) we only need to consider the completely 

symmetric case (lix,)^, and the completely antisym- 
metric case (Tx,£,)rM. 

We can further decompose the remaining subspaces by 
considering the trace of a pair of upper and lower indices. 
If such a trace vanishes, it also vanishes after applying 
the transformation of Eq. (|A3|) . and therefore such a con- 
straint defines an invariant subspace. For the completely 
symmetric case we distinguish between 



(i) 


(rp \{hj} 


= 


_v (rp \{h3} 


= 0, 


(A7) 


(ii) 








= 0, 


(A8) 


(iii) 






A (T LL )^ 




(A9) 



where we sum over repeated indices. The subspace (i) has 
nine constraints and is thus 27-dimensional, the subspace 
(ii) is orthogonal to (i) and has one constraint and is 
thus 8-dimensional, and the remaining subspace (iii) is 
1-dimensional. Therefore the completely symmetric case 
can be decomposed as 



27 < 



1 . 



For the completely antisymmetric case the analog defini- 
tion of subspaces leads to a zero-dimensional space (i) , a 
8-dimcnsional space (ii), and a 1-dimensional space (iii), 
i.e., the completely antisymmetric case can be decom- 
posed as 

0©8©1. 

This completes the classification of the left-left oper- 
ators of the form given in Eq. (|A1|) according to rep- 
resentations of SU(3)l. From now on we restrict the 
discussion to Qll with either (I) exactly one s field or 
(II) exactly one u or d field. In the case (I) only the ele- 
ments (Tll)^i and (Tll)^ with j, k, I = 1, 2 are nonzero. 
Therefore they live in a 8-dimensional representation of 
SU(2) isospin, and we can classify them according to ir- 
reducible isospin representations in the following. In the 
completely symmetric case (which is especially symmet- 
ric in k I) they live in a 6-dimcnsional isospin repre- 
sentation, and we distinguish between 



Qp7,i),i/2 = ( s d) v _ A (uu) v _ A _|_ ( su ) v _ A (ud) v _A 



where we sum over j — 1,2. The subspace (la) has two 
constraints and is thus 4-dimensional (I = 3/2), the sub- 
space (lb) is orthogonal to (la) and is thus 2-dimensional 
(I = 1/2). In the completely antisymmetric case Tll 
lives in a 2-dimensional isospin representation. The re- 
spective subspace (la) has two constraints and is thus 
zero-dimensional, the respective subspace (lb) is orthog- 
onal to (la) and is thus two-dimensional (7=1/2). In 
the case (II) the tensor transforms in the fundamental 
isospin representation (7 = 1/2). 

In the remainder of this subsection we give a list of left- 
left operators with AS = 1 transforming in irreducible 
representations of SU(3)l and isospin in which the oper- 
ator 

(sd) V - A {uu) v ^ A (A12) 
enters. The completely symmetric operators are given by 

q^27,1) : 3/2 = (gd) V ^A{uu) V -A + {su)v-A(ud) V _ A 

- (sd) v ^ A (dd) V -A , 

d)v-A{uu) V -A + 

+ 2(sd) V - A (dd) V - A - 3(sd) V - A (ss)v-A , 
q(8,1),1/2 = fsd) v _ A (uu) V -A + {su) V - A (ud) V - A 
+ 2(sd) V - A (dd)v- A 

+ 2{sd) V - A {ss)v-A, (A13) 

where (L, R) indicates the representation of SU(3)l <8> 
SU(3)r, and 1/2 (3/2) indicates the isospin. The com- 
pletely antisymmetric operator is given by 

q(8,1),1/2 = (sd) v _ A (uu)v-A - {su)v-A{ud) V - A ■ 

(A14) 



2. Left-right operators 

The left-right operators Q$ , . . . , Qg can be written as 

Qlr = (T LR )%(ql ® g*)(& ® q' R ) (A15) 

with right-handed quark fields q R . The individual quark 
fields transform in the fundamental representation of 
SU(S) L/R , i.e., under V L e SU(3) L and V R G SU(3) fl 
we find 

Q'lr = (vl) ia (V^ jb (T LR ) a c b d (V L ) ck (V R ) dl 

x (<fL®q k L)(<i R ®q R ) 

= (TL R )%(qi®qlW R ®q'k) (A16) 

with 



(la) 



(rp \{3,j} 



= o, 



(A10) {n R )l^{vl)UVl)ATL R )ti(VL)ck(V R ) dl . (A17) 
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This transformation corresponds to the 9-dimensional 
3 (X> 3 representations of SU(3)x, and SU(3)^. These 9- 
dimensional representations can be decomposed in a 8- 
dimensional subspace with vanishing trace of left (right) 
indices and in a 1-dimensional subspace with nonzero 
trace, i.e., 



3®3=8©1 



(A18) 



This completes the classification of the left-right opera- 
tors of the form given in Eq. (|A15[) according to repre- 
sentations of SU(3)l (g) SU(3)fl. From now on we restrict 
the discussion to Qlr with (Tlr) % m ^ only for (I) i = 3 
and j,k,l = 1,2 or (II) i = j = I = 3, k = 1, 2. The 
classification according to irreducible representations of 
isospin now follows from the analog discussion of left-left 
operators. 

In the remainder of this subsection we give a list of left- 
right operators with AS* = 1 transforming in irreducible 
representations of SU(3) l ® SU(3)i,; and isospin in which 
the operator 



(sd) V - A (uu) v+A 



(A19) 



enters. It is straightforward to see that Q§ and Qq trans- 
form in the (8, 1) representation of SU(3)l <8> SU(3)# and 
in isospin representation / = 1/2, i.e., we define 



n (8,l),l/2 _ n 



(A20) 



Operators symmetric under k o I are given by 

q^8,8),3/2 = (sd) v _ A (uu)y +j4 + (su) V ^ A (ud) V+A 

- (sd) V -A(dd)v+A , 

q^8,8),1/2 = (sd)y_ A (uu) V+A + {su) V - A {ud)v+A 

+ 2(sd) v ^ A {dd) v+A 

-3(sd)v-A(ss) v+ A- (A21) 

One can also construct an operator that is antisymmetric 
under k o I for k, I — 1,2 with well-defined transforma- 
tion under flavor and isospin: 

(3(8,8), 1/2 _ {gQy_ A (jiu)v +A - {su) V -A{ud)v+A 

-{sd)v- A [ss) v +A- (A22) 



lp(8,l), 1/2 _,_ \ n ( 2 



Q 9 = _±Q^l).V2 + l Q (27,l),l/2 

l n (27,l),3/2 l n (8,l),l/2 
-r 2 A 

for the color-diagonal operators and 

Q 2 = + 1 Q (27 ,l),l/2 

l n (27,l),3/2 _ l n (8,l),l/2 

3^s 2 ' 

1 „f8.1),l/2 l n (8,l),l/2 

2^A 



(A24) 



(A25) 



'4 - 



(A26) 
(A27) 

(A28) 



n _ 1 n (8,l),l/2 1 n (27,l),l/2 

^ 10 ~ - io^ 5 10^ s 

l n (27,l),3/2 _ l n (8,l),l/2 

2 s 2 A 

for the color-mixed operators with 

04= ^2 {s a qa)v-A{qbdb)v-A , 



= 2 51 e </( SQ 9a)y-A(gfcd6)y-A • (A29) 

g— u,d : s 

Operators with spinor structure (V — A)-(V + A) are 
already in explicit (8,1) and (8,8) representations, i.e., 

,i).V2 



^5,6 - ^5,6 

g7= l Q (8,8),3/2 + l Q ^8), 1 /2 : 



(A30) 
(A31) 



1 



y(8,8),3/2 1^/(8, 8), 1/2 



2^5 2~Q'a' '' ' (A32) 



with 



g,^8,8),3/2 _ ( Sa d b ) v _ A ( Ub U a )v+A 



+ (s a Ub)v-A(ubd a )v+A 

- (s a db)v-A(dbd a )v+A, (A33) 

g/(8,8),l/2 _ (5 a rf fc ) v _ j4 (i2 bUa ) y+j4 

- (SaW&V-A(M&da)v+A 

- (s Q db)\/- J 4(shSa)y+ J 4 ■ (A34) 
Appendix B: Alternative determination of AT 



3. Change of basis 

Operators of the physical basis I and II with spinor 
structure (V — A)-(V — A) can be decomposed in the 
color-diagonal operators of Eqs. (|A13J) and (|A14j) . We 
find 

Qi = -Qf 1)A/2 + -Qf ■ 1) ' 1/2 



l n (27,l),3/2 l n (8,l),l/2 



(A23) 



In this appendix we provide an alternative method to 
determine AT defined in Eq. ([54"]) relating the Wilson 
coefficients of the MS-renormalized operator basis I, II 
to the respective coefficients of the chiral basis. 

To this end we consider an arbitrary RI scheme and 
express on-shell matrix elements of the effective Hamil- 
tonian in terms of operators Q'^ 1 and first in terms of 

Wilson coefficients C'f , 
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x (q'?(aO) . 



(Bl) 



and then in terms of Wilson coefficients C 4 MS , 



G 



RI-s-MS 
ij 



(m) 



x T, 



Y?\n))- (B2) 



Since these equations hold for an arbitrary matrix ele- 
ment, we find 



MS , / " XIS / ; /j*'' 1 ' * A ' S | 1 (B3) 



where T is defined in Eq. (J58J). A comparison with 
Eq. flM} yields 



T + AT = R m ^ MS T(R 



,RI-+MSx_i 



(B4) 



Note that the left-hand side is just a conversion factor of 
Wilson coefficients in the MS scheme, and thus the right- 
hand side must also be independent of the RI scheme used 
to calculate R^™s and R ,m-+m 

We calculate AT using Eq. (|B4j) for the non- 
exceptional momentum configuration with projectors 



P, 



VV+AA,, 



(l),u 



Pi 



VV+AA,f _ VV+AA,, 

(3),u ~~ 



160iV c (^ c + 1) 



(B5) 



VV+AA, $ 



IP, 



VV+AA., 
(l),d 



2P, 



VV+AA,0\ 



(3),u 



3N r _ 



160AT C (AT C 2 - 1) 



(3P, 



(2),« 

2iV r 



2P, 



vi^+aa,^ 



16(W C (7V C 2 - l) 



VV+AA, i 
(l),u 

2N r - 



2P, 



(2),d 



VV+AA, $ 



VV+AA,*, 
(4),ti / 



(B6) 



(i),d 



2P, 



KV+AA,^\ 



(3),u 



16fW c (AT c 2 - 1) 



(3P 



VV+AA,* 
(2),u 

3N C - 2 



2P 



VV+AA,* 



(2),d 



2P, 



VV+AA,*-, 



(4), 



Pi 



160N C (N% - 1) ' 
VV-AA4 VV-AA4 



(B7) 



P, 



96 (A^ 2 - 1) 

VV-AA,* 2p VV-AA4 



(2, 



(2),<2 



P, 



96JV c (iV 2 - 1) 

1 1 -AA4 VV-AA,* 



(B8) 



P fi = 



(2),« 



(2),c2 



96 (AT 2 - 1) 



-jVV-AA, 
(l),u 



2P 



yy-AA,^ 



P 7 



96iV c (7V 2 - 1) 

W-AA,,* _ VV-AA,* 
^{l),u Ul),d 



p, 



vv 

(2), 



~ l }vv- 



p, 



AA,, 



(2),d 



Ph = 



4&N C (N? - 1) 

VV-AA,* _ VV-AA,* 
(2),u (2),d 

48(iV 2 - 1) 

VV-AA,* _ VV-AA,* 
r (l),u Ul),d 



(B9) 



(BIO) 



48A^ c (iV 2 



with 



P, -ptrcc 
i 1 4 



{Q'i 



i,j = 1,2,3,5,6,7, 



(Bll) 



(B12) 



This choice of projectors recovers the conversion matrices 
of the RI/SMOM(^, y) schemes for the non-exceptional 
momentum configuration. It makes use of projectors with 
a more complex flavor structure which allows for a defini- 
tion of this scheme without considering each irreducible 
representation of operators separately. We use the ba- 
sis of projectors {P^p AA } with i = 1, 2, 3, 4 defined in 
Eqs. The details of the wave function renormal- 

ization given in A| drop out in Eq. (|B4[) and hence we 
do not need to specify the wave function renormalization 
scheme y. The result for AT using this method agrees 
with AT given in Sec. IfflEl 
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